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Abstract 

We have searched for associated production of the lightest chargino, Wi, 
and next-to-lightest neutralino, Z2, of the Minimal Supersymmetric Stan- 
dard Model in pp collisions at -y/s = 1.8 TeV using the D0 detector at 
the Fermilab Tevatron collider. Data corresponding to an integrated lu- 
minosity of 12.5zb0.7 pb~^ were examined for events containing three iso- 
lated leptons. No evidence for W1Z2 pair production was found. Limits on 
a{WiZ2)Br{Wi lvZi)Bi{Z2 llZi) are presented. 

Submitted to the Physical Review Letters 
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Supersymmetry (SUSY) is a symmetry which relates bosons and fermions [^. Super- 
symmetric extensions of the Standard Model (SM) are attractive because they remove the 
"fine tuning" problem associated with loop corrections to the mass of the Higgs boson and 
provide a basis for gauge coupling unification at a high mass scale. One consequence of 
these models is the introduction of a SUSY partner (sparticle) for each SM state. Every 
sparticle and SM particle is assigned an internal quantum number called R-parity. If R- 
parity is conserved (as assumed in this analysis), then sparticle states are produced in pairs 
and there must be one sparticle which does not decay. This sparticle is referred to as the 
Lightest Supersymmetric Particle (LSP). The SUSY framework which introduces the fewest 
additional particles is known as the Minimal Supersymmetric Standard Model (MSSM) 0. 
If the requirement is made that SUSY be a locally invariant gauge symmetry, the result is 
a theoretical framework known as supergravity (SUGRA) 0. 

In the MSSM and minimal SUGRA there are two chargino states (Wj,i=i,2) and four 
neutralino states (.^4,4=1,4), corresponding to mixtures of the SUSY partners of the Higgs 
bosons, W and Z bosons, and the photon. (In some of the literature an alternate notation 
is used: xfi=i 2 charginos and Xij=i4 for neutralinos.) In most regions of the SUGRA 
parameter space not excluded by previous experiments, the LSP is the lightest neutralino |^ 
and thus escapes detection. The best limits to date on the masses of the Wi and Z2 
states come from the LEP experiments [Q; the current limits are > 45 GeV/c^ and 
M~ > 40 - 45 GeV/c2. 

At pp colliders charginos and neutralinos can be produced in pairs, with W1Z2 pairs 
having the largest cross section over much of the parameter space 0. Production cross 
sections (9(100-10) pb are possible at the Tevatron for Wi masses between 45 and 100 
GeV/c^ The Wi can decay into qq' or W plus an LSP, while the Z2 can decay into 

qq or II plus an LSP. The presence of neutrinos or LSP's among the decay products will 
generally lead to missing transverse energy (-^t) . The final state consisting of three leptons 
and J^rp (and little hadronic activity) has few SM backgrounds and is the subject of the 
present analysis. 
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The spectra of the transverse momenta (pt) of the final state leptons can be relatively 
soft due to the three-body decays of the Wi and Z2 involving massive non-interacting 
particles. Figure |I] shows the expected pr spectra of the final state leptons as well as 
the ]^rp distribution at the physics generator level for simulated W1Z2 31 events, with 
= 56 GeV/c^. These Monte Carlo events follow the mass relation common to many 
SUSY models: M,j. ^ ^ 2M^ . 

The data used in this analysis were obtained using the D0 detector at the Fermilab 
Tevatron pp collider operating at a center of mass energy of 1.8 TeV. The total integrated 
luminosity used in this analysis from the 1992-1993 Tevatron run was 12.5 ± 0.7 pb~^. 

The D0 detector has three major subsystems: central tracking detectors (with no cen- 
tral magnetic field), uranium-liquid argon electromagnetic and hadronic calorimeters, and 
a muon spectrometer. The detector and data acquisition system are described in detail 
elsewhere 0. The central tracking system is used to identify charged tracks in the pseudo- 
rapidity range |?7| < 3.5. The calorimeters provide full angular coverage for |?7| < 4.0, with 
transverse segmentation A?7 x A0 = 0.1 x 0.1, where is the azimuthal angle. The muon 
system consists of proportional drift chambers and magnetized iron toroids with coverage 
extending to |?7| < 3.3. 

Electrons were identified as calorimeter clusters having at least 90% of their energy 
deposition in the electromagnetic calorimeter, with one or more tracks pointing to the cluster. 
Jets were reconstructed from energy deposition in the calorimeters using a cone algorithm 
with cone size IZ = \^Ar]'^ + = 0.5. Muon tracks were reconstructed using hits in the 
muon drift chambers; their momenta were calculated from the bend of the tracks in the 
toroids. 

Combinations of single lepton and dilepton triggers were used for the four final states 
(eee, ee/i, efifi, and fififi). These triggers included: a single muon with pj, > 15 GeV/c; two 
muons with pj. > 3 GeV/c; one muon with pt^ > 5 GeV/c plus one electromagnetic cluster 
with > 7 GeV; one electromagnetic cluster with > 20 GeV; and two electromagnetic 
clusters with E^ > 10 GeV. The integrated luminosity per channel is given in Table |. 
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Events passing the trigger requirements were selected offline by requiring three or more 
reconstructed leptons (electrons or muons) having > 5 GeV or pi^ > 5 GeV/c, with 
|?7e| < 2.5 or |?7^| < 1.7. There were 2827 events in this initial data sample. Electrons and 
muons in these events were then required to pass the quality cuts described below. 

Electrons were required to have transverse and longitudinal shower profiles consistent 
with expectations based on detailed Monte Carlo studies [|lOl, to have no more than two 



tracks pointing to the calorimeter cluster, and to have an electromagnetic isolation X < 0.15, 
where X = [Etot — Eem]/ Eem, -E'tot is the total cluster energy inside a cone of radius 71 = 0.4, 
and Eeu is the electromagnetic energy inside a cone of 7^ = 0.2. For electrons with Et 
between 5 and 10 GeV, the isolation cut was relaxed to X< 0.2 to increase efflciency. 

Muons were required to have a separation from any jet of at least TZ = 0.5, to be aligned 
with minimum ionizing energy deposition in at least 50% of all calorimeter layers and in 
at least 60% of the hadronic calorimeter layers, and to have either a matching track in the 
central detectors or impact parameters in the rz (bend) and xy (non-bend) views consistent 
with the muon having been produced at the primary event vertex |T^. To reduce cosmic 
ray background, muons were required to be in time with the beam crossing and any muon 
pair having both polar and azimuthal opening angles greater than 165° was rejected. 

There were 19 events after these quality cuts. The following topological cuts were applied 
to these events. For the eee channel, events were required to have > 10 GeV, with the J^j^ 
reconstructed using only energy deposited in the calorimeters. This cut reduced background 
from Z/'-f ^ e~^e~ events with a third electron from either a photon conversion (including 
vr^ 77) into an unresolved e~^e~ pair or a jet which was reconstructed as an electron. 
Since extra material in the forward region enhances the photon conversion probability, the 
data exhibit an excess of electrons in the forward region while the signal distributions peak 
in the central region. Therefore, a cut was applied in the eee and ee/i channels to exclude 
events with more than one electron in the region Ir^l > 1.7. For the efifi and fififi channels, 
muon pairs were required to have an invariant mass greater than 5 GeV/c^, which reduced 
background from J/ip events and the combinatoric background in the reconstruction of 
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muons. Table | summarizes the effect of the cuts on each of the channels. We see no 
candidate events consistent with W1Z2 pair production and subsequent decay into trilepton 
final states. 

Detection efficiencies were determined using a combination of data and Monte Carlo 



simulations. Monte Carlo signal events were generated using ISAJET [0] and processed 



with a full simulation of the D0 detector based on the GEANT program. Seven sets of 
events were generated, with the mass of the Wi varying from 45 to 100 GeV/c^. Because of 
the correlation between the masses of the Wi, Z2, and Zi, efficiencies can be parametrized 
as a function of M{r,. These Monte Carlo events were used to determine kinematic and 
geometric acceptances only. 

Electron identification efficiencies were determined from a set of simulated single electron 
events generated in six bins between 5 and 25 GeV. These were overlaid with minimum 
bias events from collider data in order to include the effects of the underlying event and any 
noise in the calorimeter on electron isolation and shower profile. The results of these studies 
for high Et electrons were verified by analyzing a sample oi Z ^ ee events in which one 



electron was required to pass all cuts and the second electron was then used as an unbiased 
estimator for each cut. 

Similarly, muon identification efficiencies were based on Z fifi and J/ip ^ /^/^ event 
samples. These two sets provided independent estimates of efficiencies for both high and 
low Pt muons. 

Electron and muon identification efficiencies were parametrized as a function of the 
electron E^ or muon px and incorporated with the topological cuts described above to 
determine the overall analysis efficiency for each set of Monte Carlo signal events. These 
efficiencies are shown in Fig. |^ for each final state, along with a parametrized fit |T^, as a 
function of the Wi mass. 

Backgrounds were estimated from data whenever possible, supplemented with Monte 
Carlo simulations. Standard Model processes which produce three or more isolated leptons, 
such as vector boson pair production and semileptonic decays in heavy flavor production, 
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are expected to yield less than 0.1 event in any channel. Thus the primary sources of 
background are single lepton and dilepton events with one or more spurious leptons. The 
sources of spurious electrons are jet fluctuations and unresolved e^e" pairs from photon 
conversions. The probability of a jet faking an isolated muon is negligible. 

The background from fake electrons was calculated from data using dilepton events with 
one or more additional photons and/or jets. The expected number of events was determined 
by multiplying the number of events seen in data by the probability of a photon conversion 
or the rate for a jet to fake an electron [T^. The primary source of background in the ixixfi 



channel is heavy flavor {hb and cc) events with the muons produced at large angle to the 
jets. The total background for each flnal state is included in Table |. 

Based on zero candidate events, we present a 95% confidence level upper limit on the 
cross section for producing W1Z2 pairs times the branching ratio into any one of the trilep- 
ton final states. The results from the four channels were combined in the calculation of the 
limit, with the assumption that Br(eee) = Br(ee/i) = Br(eyU/i) = Br(/i/iyu). Uncertainties 
in this calculation include the uncertainty in the luminosity (5.4%) and uncertainties in the 
overall analysis detection efficiencies (between 15% and 25% of the value) due to Monte 
Carlo statistics, systematic errors in the determination of lepton identification efficiencies, 
systematic errors in the trigger efficiencies, and systematic errors arising from energy scale 
corrections. To construct this limit we used the Bayesian approach of 0], with the distribu- 
tion of systematic errors represented by a Gaussian and a fiat prior probabihty distribution 
for the signal cross section. 

In Fig. 1^ we show the resulting limit in the region above the LEP limit 0]. For compari- 
son, we show three bands of theoretical curves. Band (a) shows the ISAJET production cross 
section obtained with a wide range of input parameters, multiplied by a branching ratio of 
|. The value of | for a single trilepton channel is obtained when the Wi and Z2 decay purely 
leptonically and lepton universality is applied. Branching ratios of this order are predicted 
in models with very light sleptons, as for example the model of Ref . . Bands (b) and (c) 
show the ax Br values from ISAJET obtained with the following SUGRA input parameters: 
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mo = [200,900] GeV/c^, mi = [50, 120] GeV/c^ A = and the sign of // negative. Band 
(b) is for tan 13 — 2 and band (c) for tan /3 = 4. 

In conclusion, we have searched for the associated production of chargino and neutralino 
pairs by looking for the reaction pp W1Z2 31 + X. We see no evidence for W1Z2 pro- 
duction in 12.5 pb"^ of data. This leads to upper limits on cr (W1Z2) Br ii/Zi)Br(Z2 
llZi) ranging from 3.1 pb for M^^ = 45 GeV/c^ to 0.6 pb for M^^ = 100 GeV/c^. 
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Channel 


eee 


ee/x 


e/x/i 


uuu 


JCdt (pb-i) 


12.5 


12.5 


12.2 


10.8 






Events Remaining 








By Analysis Channel 






13 


42 


297 


2475 


With Quahty Cuts 


5 


2 


5 


7 


Ne forward < 2 


4 





N/A 


N/A 


$;t > 10 GeV 





N/A 


N/A 


N/A 


M^^ > 5 GeV/c2 


N/A 


N/A 








Candidates 














Background 


0.8 ±0.5 


0.8 ±0.4 


0.6 ±0.2 


0.1 ±0.1 



TABLE I. Analysis cuts for each of the search channels, showing the number of events left 
after a cut has been applied. No candidates are seen in any of the four channels. The predicted 
background per channel is also shown. 
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